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Abstract 


The  problems  of  pressure  shock  and  longitudinal  impact  on 
a  semi-infinite  circular  tube  are  studied.  The  considerations 
are  performed  within  the  framework  of  the  linear  theory  of 
elasticity.  The  problems  were  reduced  to  the  system  of  Bessel 
differential  equations  by  using  double  integral  transforms.  The 
solutions  are  obtained  in  the  terms  of  Fourier  and  Laplace 
integral  which  could  be  evaluated  by  numerical  methods. 
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Introduction 

Problems  of  wave  propagation  in  waveguides  in  view  of 
their  importance  in  numerous  applications  have  received 
considerable  attention  in  the  literature  devoted  to 
elastodynamics .  The  bulk  of  the  study  of  elastic  waves  is 
concerned  with  rods  and  plates.  The  cylindrical  rod  of  circular 
cross  section  is  a  classical  example  of  a  medium  involving  the 
complexities  of  dispersive  propagation.  Pochhammer  [1]  was  the 
first  who  derived  frequency  equations  for  the  cases  of 
compressional ,  flexural  and  torsional  harmonic  waves  traveling 
in  an  infinite  elastic  rod.  The  frequency  equations  are 
connected  with  a  stress  free  lateral  surface  of  the  rod  and 
they  constitute  a  basis  for  investigations  of  both  steady 
vibration  and  transient  wave  propagation. 

This  report  deals  with  two  dynamic  problems  of  a 
se.mi-inf inite  circular  elastic  tube.  The  first  problem  is 
related  to  mixed  pressure  conditions,  the  second  one  is  devoted 
to  the  longitudinal  impact.  The  analysis  is  performed  within 
the  framework  of  the  linear  theory  of  elasticity.  The 
considered  problems  belong  to  the  class  involving  mixed  edge 
conditions . 

The  method  of  solution  of  the  first  problem  is  based  on 
the  double  integral  transforms  (Fourier  and  Laplace  transform  ) 
and  it  leads  to  solutions  given  in  the  terms  of  Fourier  and 
Laplace  integrals  ,  which  can  be  calculated  by  using  analytical 
approximations  or  numerical  analysis. 

A  similar  mixed-pressure  problems  of  infinite  circular 
rods  were  solved  by  Folk  et  al  [2]  and  is  also  presented  in  the 
Achenbach  book  [3].  By  using  integral  transforms  and  employing 
the  analysis  of  Pochhammer's  frequency  equations  to  the 
procedure  of  an  inversion  of  Laplace  transforms,  the  solution 
of  the  wave  problem  was  expressed  in  the  terms  of  Fourier 
integrals.  The  Fourier  integrals  could  not  be  evaluated 
rigorously  by  analytical  means,  however,  a  satisfactory 
approximation  was  obtained  by  the  method  of  stationary  phase. 
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The  second  problem  considered  in  this  report  deals  with 
the  impact  of  a  semi- inf ini te  circular  elastic  tube.  The 
solution  is  given  under  assumptions  that  the  end  of  tube  has  no 
tangential  stresses  but  is  suddenly  set  into  motion  with  a 
constant  velocity  and  the  lateral  surfaces  are  free  of 
pressures.  The  analogous  problem  for  a  semi-infinite  bar  was 
solved  by  Skalak  [11].  His  method  was  to  split  the  problem  into 
two  parts.  In  the  first  part  the  collision  of  two  half-spaces 
was  considered.  In  the  second  part  he  investigated  an  infinite 
bar  with  an  equal  and  opposite  stresses  travelling  along  the 
lateral  surface.  Using  the  principle  of  superposition  a  stress 
free  lateral  surface  was  ensured.  In  Section  2  Skalak' s  results 
are  adopted  to  the  problem  of  circular  tube.  The  solution  of 
the  impact  problem  is  presented  in  the  form  of  Fourier 
exponential  and  Laplace  integrals,  which  can  be  calculated  by 
using  analytical  approximations  or  numerical  analysis. 

Transient  wave  problems  in  elastic  waveguides  were 
considered  by  many  authors.  In  the  survey  by  .Miklowitz  [4]  one 
can  find  a  coverage  of  general  literature  up  to  1964.  Solutions 
of  many  wave  problems  are  also  presented  in  Miklowitz  book  [5]. 
Longitudinal  harmonic  motions  of  a  cylindrical  tube  on  the 
basis  of  the  three-dimensional  theory  of  elasticity  were 
considered  by  Gosh  [6],  where  the  frequency  equation  are 
derived.  The  discussion  of  the  characteristic  equations  for 
this  problem  and  the  phase  velocity  (for  the  case  of  Poisson's 
ratio  r*  -  0.3)  are  presented  in  [7].  Some  numerical  results  for 
the  longitudinal  impact  problem  of  a  circular  tube  obtained  by 
using  the  numerical  integration  along  the  bicharacteristics, 
are  given  in  [8]  . 

Understanding  of  dispersion  behaviour  in  bars  and  tubes 
are  of  great  importance  in  all  experimental  techniques  used  to 
test  mechanical  characteristics  of  materials  at  high  loading  or 
strain  rates,  for  example  [12-14].  One  prominent  example  is  the 
Split  Hopkinson  Bar  technique  (Kolsky  apparatus),  where  a  thin 
wafer  specimen  is  inserted  between  two  instrumented  Hopkinson 
bars  [12].  A  more  exact  analysis  of  wave  dispersion  in  the  SHB 
enables  for  a  more  exact  interpretation  of  oscillograms,  and 
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consequently,  for  a  more  exact  determination  of  a  specimen 
response  at  different  loading  conditions  [15,16]. 

Since  it  is  expected  that  the  lateral  inertia  causing  the 
Pochhammer-Chree  vibrations  are  substantially  reduced  in  tubes, 
application  of  tubes  instead  of  bars  as  measuring  devices  in 
wave  transmission  seems  to  be  reasonable.  So  far,  some 
applications  of  tubes  for  material  testing  at  high  strain  rates 
were  already  proposed,  [17,18].  However,  in  those  applications 
a  simple  one-dimensional  wave  appro:  imat ion  was  commonly  used. 
The  main  target  of  this  study  is  to  make  analysis  of  the 
wave-guide  tubes  more  precise,  specially  for  short  pulses. 

1 ■ The  mixed  problem  of  pressure  shock 

1.  1  Formulat ion  of  the  problem 

A  semi-infinite  cylindrical  elastic  tube  of  inner  radius 
a  and  outer  radius  b  ,  which  are  referred  to  cylindrical 
coordinates  {r,&,z)  where  the  z-axis  coinciding  with  the  axis 
of  the  tube,  Fig.i,  is  considered.  Let  X,  p  be  Lame  constants 
and  p  be  the  mass  density  of  the  tube  material.  The  problem  is 
restricted  to  the  axisymmetric  case  of  motion,  thus  the 
displacement  vector  u  is  independent  on  the  variable  &  and 

u(r,z,t)  =  (u(r,z,t);  0;u(r,z,t)  )  (1-1) 


linear  elasticity  take  the  form 


dr  dz 


r  z 

where  t  denotes  the  time. 
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The  stress  components  are  expressed  by  the  displacements 
u  as  follows 
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Consider  the  boundary  conditions 
(  r,  z  =  0,t  )  =  -  f  (t) , 


(  r,z  =  0,t  )  =  0  ,  for  r  e  <a,b>  ,  t  *=  (0,od)  (1.4) 

and 


a  (r=a,z,t)  =  o  (r=b,z,t)  =  0  , 

rr  rr 

a  (r  =  a,z,t)  =  a  (r  =  b,z,t)  =  0,  for  z  e  (O.co)  ,  t 

rz  rz 


(1.5) 
(0,oo)  , 


where  f(.)  is  a  given  function,  provided  that  f(t)  h  0  for  t<0. 
A  futher  condition  is  that  the  displacements  and  the  stresses 
vanish  at  infinity  ,  z  — >  <»  .Assuming  that  the  semi - inf inite 
circular  tube  is  at  rest  prior  to  time  instant  t  =  0,  equations 
(1.4)  and  (1.5)  are  supplemented  by  the  initial  conditions 
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du 


u^  (  r ,  z ,  t=0 ) 

=  — ’’  (  r ,  z  ,  t=0  )  =  0  , 
dt 

(1.6) 

u^(r, z, t=0) 

=  — — ^(r,z,t  =  0)  =  0  , f or  r  e  <a,b>  ,z 
dt 

€{  0,00)  . 

By  introducing  the  displacement  potentials 

as  follows 
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(1.7) 

the  equations  of  motion  (1.1)  are  separated  into 
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(1.10) 
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Substituting  equations  (1.7)  into  (1.3)  the  stress  components 
may  be  written  in  the  terms  of  displacement  potentials  in  the 
form 
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From  boundary  condition  (1.4)^  and  equation  (1.7)^  it 
that 

d<t> 


dr. 


dz 


,  for  z  =  0  ,  re  <a,b>  ,  t  «  (0,00). 


f ol lows 


(1.13 


By  using  (1.13)  ,  (1.12)^  and  boundary  condition  (1.4)^  one 
obtains 

d^4> 

p  -  (r,z=0,t)  =  -  f(t)  ,  for  r  e  <a,b>,  t  e  (0,®)  .  (1.14) 


Because  the  potential  <P  is  independent  on  r  at  z  =  0  (see 
Eq.(1.14)  ),  then  d4>/dx  =  0  for  z  =  0  ,  and  from  (1.13)  it 
follows  that 


—  (r,z=0,t)  =  0  for  r  e  <a,b>  ,  t  e  (0,®)  .  (1.15) 

dz 

Equations  (1.14)  and  (1.15)  constitute  boundary  conditions  at 

^  ^  ^  ^  ^  4  ^  ^  ^  (— %  ^  2  ^  *— t  ^  ..  A-  •% 

Z  —  u  JL  Ox.  t-itc;  ^  X  X  u  c.t:iixi.xCLx&-  jiiAC  m  xl  m  cx  ^ 

(1.6)  together  with  (1.7)  imply 
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d<p 

(p{T  ,z,Q)  =  —  (  r,  z  ,  0 )  =  0  , 

dt 

V'(r,z,0)  =  —  (r,z,0)  =  0  .  (1.16) 

dt 

Moreover,  the  potentials  4>  and  has  to  vanish  at  infinity 

,  Z  — >00  . 


1.2  Solution  by  integral  trans/orns 


Let  us  denote  by  .according  to  [9]: 

00 
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®  o 


00 
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and 
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£  {  — --  ;  t->p  }  =  f(r,z,p)  -  p  iim  f(r,z,t) 
at  t->o 

af 

lim  -  ( r, z, t )  (1.19) 

t->o''  at 

provided  ,  that  |f(r,z,t)|  <  K(r,z)  exp(?"^t)  ,  where  K(r,z)  is 
an  arbitrary  function  and  Re  p  =  ?-  >  ^^  • 

Applying  Fourier  sine  transform  with  respect  to  z  to 

eauation  (1.8)  it  follows  that 
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a^<p^  d<p^  1 

— +  - - -  ?  0(r,z  =  O,t)  =  (1.20) 

ar  ^  ar  C  at 


and  using  Laplace  transform  to  (1.20) 
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(1.21) 

The  boundary  condition  (1.14)  after 

yields 

Laplace 

transformation 

-C  f ( P ) 

<P  (r,z-0,p)  =  -  , 

(1.22'' 

pp 

where 

f(p)  =  -i?  {  f(t)  ;  t  ->p  }  . 


(1-23) 


Substituting  (1.22)  into  (1.21)  the  following  Bessel  equation 
is  found 


a^^^(r,f ,p) 


ar 


1  a0^(r,?,p) 

^  ar 


1- 


(p) 


(1.24) 


PP 


where 
a  =  ^  + 


(1.25) 
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Applying  Fjurier  cosine  transform  with  respect  to  z  to  equation 
(1.8)^  one  finds 

1  dv  1  1 

— ^  - _  _  (r,z  =  0,t)  -  ■=  _  _  .  (1.26) 

at'  ^  dr  a-L  r  c  ax 

2 

However,  by  means  of  boundary  condition  (1.15)  the  equation 

(1.26)  is  reduced  to  (after  Laplace  transform): 

^2  — c  .  -~c  , 

a  14)  \  a^p  1  _ 

+  - -  -  {(5^  +  -  )  =  0  ,  (1.27) 

ax  ^  dr  r 


where 


=  k""  * 


(1.28) 


The  general  solutions  of  equations  (1.24)  and  (1.27)  may  be 
written  in  the  form,  [5]:^ 


{r,?,p)  =  AJ^(a  r)  +  A^K^(a  r)  -  . f(p)  , 

p  p  a 


p'^(r,^,p)  =  B^I^(/7  r)  +  B^K^(/?  r)  ,  (1.29) 

where  (  -  )  /  ^  ^  modified  Bessel 

functions  of  the  first  and  the  second  kinds,  respectively.  Here 
,  A^,  are  unknown  functions  of  the  variables  p  and 

they  have  to  be  determined  by  the  use  of  boundary  conditions  on 
the  lateral  surfaces  of  the  tube  given  in  (1.5).  In  order  to  do 
this  the  equations  (1.12)  are  transfomed  ,  thus 

a^p^  dhp^  p^  __ 

=  2r^  (  <p^  .  (1.30) 

dr  <?r  c 


^It  is  emphasized  that  for  the  analogy  problem  of  a  rod,  the 
following  conditions  should  be  assumed  :A  =B  3  0  to  ensure 
finite  potentials,  displacements  and  stresses  at  r=0. 
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From  (1.15)  ,  (1.13)  and  (1.31)  it  is  obtained 


d0“ 

— C  ,  K  «2— C  ,  2  C 

O-  -  2/u  (  ?  -  -  +  K  ^  )  +  p  p  V' 

di¬ 


ll  .  3  3  ) 


Employing  the  following  relations,  according  to  [10]: 
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and  using  equation  (1.29)  it  is  found 
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(1.34) 


Substituting  equations  (1.29)  and  (1.34)  into  (1.30)  and  (1.32) 
one  finds 
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(r,?,p)  =  2p  [  a  IJcxr)  -  o.  K^(ar)  +  B^{?  + 


^  )  I,(^r)  .  B^(?  .  )  K  (^r)  ] 


(r,?,p)  =  A  [  {2u  +  — _-  )  I-(«r)  -  2(j  I,  («r)  ]  + 

r  r  1  ^  u  L  * 

C 

\  2  1 

Xp 

A^  [  (2u  +  — -  )  K^(o<r)  +  ( ar )  ]  + 

c 

1 

2:j  K  [  ft  l^Wr)  -  ~  IJ^r)  ]  +  2iu  ^  [~ft  K^iftt)  - 


i  K  {ftv)  ]  - -  f  (P)  . 

^  (X  +2p)c« 


(1.35) 


Let  us  introduce  the  following  notations: 
e^  {  r)  ^  K  c<  (ar )  , 


e2(r) 


-  ?  c*  K^(ar)  , 


e3(r)  =  (?^  +  -2—  )  I^(/?r)  , 


e_^(r) 


..2 

If  + 


K^(^r)  , 


(1.36) 


d^(r)  =  (2aj  o<^  +  — ■-  )  -  2aj  p  Ij(cxr) 

c 

1 


d^Cr)  =  (2^  a*  +  — -  j  K^(«r)  +  2a/  p  K^(o«r) 

c 

1 


d^Cr)  =  2a/  f  {ft  l^inr)  -  p  ]  . 


i^(r)  s  -2p  f  [  /?  K^(r?r)+  p  \iftT)  1  . 


The  boundary  conditions  (1.5)  together  with  equations  (1.35) 
and  (1.36)  yield  the  system  of  four  linear  algebraic  equations. 


11 


for  unknowns  A  ,  A  , 
1  2 

A^  d^(a)  +  A^  03^(3) 

A^  d^(b)  +  A^  d^(b)  • 

A^  e^{a)  +  A^  e^ia)  - 


B  ,  B  ,  in  the  form: 
1  2 

^  ®2 

B^  d^(b)  +  d^(b) 

+  ^2 


^  ?  f{p) 

(X  +  2a^)  oi^ 

^  ?  f(p) 

(X  +  2a^)  ol^ 

0  ,  (1-37) 


B  e  (b)  +  B  e  (  d) 

1  3^  '  Z  4  ' 


0  . 


The  solution  of  equations  (1.37)  can  be  written  in  the  form 


where 


(1.38) 


dj(a) 

d^(a) 

dgCa) 

d^(a) 

d,(b) 

d^Cb) 

d3(b) 

d^(b) 

A  =  det 

e^{a) 

e^Ca) 

63(3) 

e^{a) 

e^(b) 

e2(b) 

63  (b) 

e^(b) 

i  »  ^ 

3 

4  ' 

i 

^  K  f(p) 

(X  +  2m  )at^ 

d^(a) 

d3(a) 

d^(a) 

=  det 

^  ?  fTp) 

(\+2p)a^ 

d^(b) 

daCb) 

d^(b) 

0 

ej,(a) 

63(3) 

e_^(a) 

0 

e^(b) 

SgCb) 

e_^(b) 

12 


X  ?  f{p) 
(X  +  2aj  )c<^ 
X  ?  f(p) 


dgCa) 


(  a ) 

4 


W  =  det 

2 

djb) 

(X+2p )a^ 

cl3(b) 

d^(b) 

e^(a) 

0 

e3(a) 

e_^(a) 

ejb) 

0 

e3(b) 

e^(b/ 

(1.39) 


W  =  det 

3 


X  ?  f(p) 


d^(a) 

d^(a) 

(X  +  2aj  )c»^ 

d^(a) 

d^(b) 

d^Cb) 

X  ?  f(p) 

(\+2p)a^ 

d^(b) 

e^{a) 

e^Ca) 

0 

e  (a) 

4 

e^{b) 

e2(b) 

0 

e  (b) 

4 

W  =  det 

4 


d^  (a) 

d^(a) 

d3(a) 

d^  (b) 

d2(b) 

d3(b) 

(a) 

Sj,  (a) 

63(3) 

(b) 

e2(b) 

63  (b) 

X  t  f{p) 
(\  +  2/l/  )c^ 
X  ?  f(p) 
(X  +  2/J  )o^ 
0 


Substituting  the  solution  given  by  (1.38)  into  (1.29)  one  has 

W,  ?  f(p) 


0  (r,? ,p) 


W  W 


A  ^ 

A  o 


Z  2 
p  p  a 


V/'^(r,?,p)  =  +  a”*  - 


(1.40) 
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From  equations  (1.7)  and  (1.40)  it  follows  that 


u  ( r , p) 


Wi 

=  -  a  I^(ar)  -  ^  o.  K^(o<r)  +  ^  ^ 

W 

?  K^((?r)  , 


(1.41) 


-c  '^1  ^ 

u^(r,?,p)  =  J  ^  A  ^  K^(«r)  ^  0 


W 


n  Kint)  + 


?  f(p) 


p  p 


Substituting  (1.38)  into  (1.35)  one  obtains 

-c  2a. 

^2,(r,?,p)  =  ^  a  I^(ar)  -  ?  a  (ar  )  + 


PP 


PP 


27;  )  I,{(5r)  f  (  ?  +  2-  )  K^(A3r)  ]  . 


W 


=  2:'  [  (2P  )  Io(o.r)  -  2u  ?  IJar)  ] 


^2  .  3 


[  (2/j  a  +  —  )  K  (ar)  +  2p  -  K  (ar)  ]  + 

o  r*  1 


w 


1  w 

2P  ?  [  /I  lo(^r)  -  -  I^(/7r)  ]  -  2^  !^  [  ft  K^(^r)  f 


7  K^(A?r)  ]  - 


^  ?  f(p) 

(X  +  2A.)a^ 


(1.42) 


Using  equations  (1.12)^  and  (1.40)  after  some  manipulations  one 
has  ,, 

— s  W  W 

o'^^(r,?,p)  =  (  pf?  -  2p  «^)  I^(«r)  +  -^  (  p 


2p  a  )K^(ar) 


W 

■g  2p  K  0  Io{^r)  + 


W 


g  2p  K  ft  K^iPr)  -  ~  f(p)  . 

a 


(1.43) 
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Equations  (1.41)  ,  (1.42)  and  (1.43)  represent  the  solution 
(called  the  formal  one)  of  the  considered  problem  expressed  in 
the  terms  of  Fourier  sine-Laplace  or  Fourier  cos  ine-Lapiacfj 
transforms.  It  should  be  emphasized  that  'i,  a  ,  /li  ,  w 
v.  =  l,2,3,4,  (see  Eqs .  (1.25),  (  1.28)  ,  (1.36)  and  (1.39)  )  are 

functions  of  transform  parameters  ?  and  p.  Using  the  inversion 
theorems  of  the  Laplace  and  Fourier  transforms  the  solution  can 
be  written  in  the  form 


u  (  r ,  z  ,  t ) 


u^  (  r ,  z  ,  t ) 


{r,z,t)  = 


1 

2  . 


■S’  (  r ,  z ,  t )  =  — ~ 

rr  '  '  '  2  . 


CD 

j^  +  ioo 

s 

S 

u^( 

0 

ioo 

r 

00 

^•  +  ioo 

S 

S 

U^i 

0 

r-iao 

X 

oo 

y  +  ioo 

f 

S 

X  r 

0 

r-ico 

CO 

iCB 

S 

S 

--'S 

a 

0 

y-ioa 

r  r 

00 

y  +  ioo 

/ 

/ 

--s 

<7 

0 

j'-  ioo 

X  X 

if 

the 

.Pt 


,Pt 


sPt 


(1.44) 


(1.4)^,  is  yet  undefined.  The  possibility  of  the  following 
cases  of  impulses  can  be  mentioned: 


Case  i  In/ini t&Ly  short  impulse 


Let  function  f(.)  (see  Eq.(1.4)^  )  be  given 

f(t)  =  -f^  6(t)  (1.45) 

where  f  ^ ,  f^  >0  ,  is  a  constant  and  <5(.)  is  the  Dirac 

function,  then 

f(p)  =  -  f^  .  (1.46) 
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Case  2  Ir\f  ini  tely  long  imfrulsG 


It  is  assumed 


f{t) 


f  H(  t 
o 


(1.47) 


where  H{ . )  is  the  Heaviside  step  function,  then 


f{p)  =  f^/p 


(1.48) 


_Cas^  3  Impulse  o/  finite  lenght 


It  is  assumed 


f(t)  =  -f^  H(t)  H(t  -  t 


(1.49) 


where  t,  t  >  0,  is  a  duration  of  the  impulse.  In  this  case  one 
has 


f(P)  =  -  p°  {  1  - 


(1.50) 


2 ■ The  problem  of  longitudinal  impact 

The  problem  presented  above  is  related  to  the  mixed 
conditions  of  shock  pressure  at  the  end  of  the  semi - inf ini te 
circular  elastic  tube.  Consider  now  the  problem  shown  in  Fig. 2. 
Two  semi-infinite  circular  tubes,  moving  in  opposite  directions 
with  speed  v  , are  assumed  to  contact  at  the  instant  t=0  and  at 
the  plane  z=0.  The  lateral  surfaces  of  tubes  are  free  of 
loads.  The  analogical  problem  for  two  semi-infinite  circular 
bars  was  solved  by  Skalak  [11],  and  such  approach  and  also  some 
results  of  [11]  will  be  used  further  on. 

The  solution  of  the  problem  considered  is  derived  by 
application  of  the  superposition  principle  for  two  separate 
problems.  The  first  of  them  is  limited  to  two  elastic 
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semi-spaces,  moving  in  opposite  directions  with  speed  v.  They 
are  assumed  to  contact  at  the  instant  t=0  and  at  the  plane  z~0. 
After  time  t=0  the  semi-spaces  are  assumed  to  behave  as  a 
single,  solid  space,  like  the  tubes  considered  above,  which  are 
assumed  to  behave  for  t  >  0  as  a  single,  infinite  tube.  The 
solution  of  this  problem  is  given  in  [11]  in  the  form: 


u^  (  r ,  z  ,  t )  =  0  , 


u^  (  r ,  z ,  t ) 


+  v  t 
-  V  z/c^ 


for  -  z  >  c  t^ , 
for  I z I  <  c^  t 


(2.1) 


The  displacements  assumed  above  induce  the  following  stresses, 
when  using  equations  (1.3)  and  (2.1): 

O'  {  r ,  z  ,  t )  =  0  , 

zr 


(  r  ,  z  ,  t ) 

^„{r,z,t) 


-{X  +  2p) 

(  “ 

\  -  X  V  /c 


for 

V  /  c  for 

for 

for 


Izl  >  c^t 
|z|  <  c^t 

|z|  >  c^t  . 
\z\  <  c  t  . 


(2.2) 


It  is  seen  from  equations  (2.2)  that  the  solution  of  the  first 
proble.m  cause  non-zero  radial  tractions  on  the  lateral  surfaces 
of  the  circular  tube.  Applying  the  principle  of  superposition 
to  ensure  a  stress  free  lateral  surfaces  of  the  infinte  tube 
the  second  problem  is  considered  ,  which  is  defined  by 
equations  of  motion  (1.2),  the  boundary  conditions 


cx^^(r=a,z,t) 

^rr(r=a,z,t) 

for  z  €  (-00 


=  r=b,  z,  t )  =  0  , 
=  <",,(r  =  b,z,t)  =(  ^ 

,  +00)  ,  t  e  ( 0 ,  +00) 


v/c. 


for  |z|  >c^t 
for  |z|  <c^t 


(2.3) 


and  the  initial  conditions  given  in  (1.6)  for  z  e  (-oo.+oo).  A 
futher  condition  is  that  the  displacements  and  the  stresses 
vanish  when  |z|  ~>  oo  . 
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The  solution  of  the  problem  is  found  by  using 

displacement  potentials  and  V'  introduced  in  (1.7)  and  Fourier 
exponential  transforms 

fMr,?,t)  H  ^{f{r,z,t);  z->?  }  =  S  f{r,z,t)  e  ^  dz  (2.4) 

-00 


as  well  as  Laplace  transform  defined  by  (1.17)^  with  respect  to 
t  .  Using  (2,4)  one  has 


af 

Jz  ' 

z->?  }  =  i  ?  f^(r,f,t)  , 

a^f 

az^ 

;z->?  }  -  f^(r,f,t) 

(2.5 

^f ( r , z , t ) 

p  rovided,  that  lira  f(r,z,t)  =  lim  a  '  “  ^  • 

I  z  I  ->oo  I  z  I  —>00  ^ 


Applying  Fourier  exponential  transform  and  Laplace  transform  to 
equations  (1.8)  one  obtains 

a^^^(r,f,p)  .  a0^(r,?,p)  2  _,p. 

- . . , -  +  -  -  -  a  <;!>(r,?,p)=0, 

ar  ^  ar 

a^v/^(r,?,p)  .  av^{r,?,p)  , 

- - -  +  i  - - - -  i  )  v^‘-(r,f,p)  =  0  (2.6) 

ar  ar  r"' 

2  2 

where  a  is  given  in  (1.25)  and  ft  is  determined  by  (1.28)  . 

The  general  solution  of  equations  (2.6)  ^an  be  written 

2 

in  the  form: 

0^(r,?,p)  =  ^  ^2  ' 

<;/^(r,^,p)  =  I^(ar)  +  K^{ar)  ,  (2.7) 

where  C  ,  C  ,  D  ,  D  are  unknown  functions  of  the  variables  K, 

1  2  1  2 

p  and  they  have  to  be  determined  by  using  boundary  conditions 

(2.3).  Applying  Fourier  and  Laplace  transforms  to  equations 

It  has  emphasized  that  for  the  analogic  problem  of  a  rod  the 

conditions  C  =D  =  0  sholud  be  assumed  to  ensure  finite 

2  2 

potentials,  displacements  and  stresses  at  r=0  . 


(1.12)  one  has  (by  using  (2.5)  and  (1.19)  )  : 

P 

O'  ( r ,  f  ,  p )  =  2g  (  -  i  ?  -  )  +  'P 

ar  ar  c 

1 

-F  2  -F 

ol  (v,K.9)  =  2p  [  i  ?  - +  (e  +  -^^-  )  V'  ]  .  (2.8) 

ar 

Substituting  (2.7)  into  (2.8)  and  employing  (1.33)  it  follows 
that 

-F  2 

[  (2p  a  +  )  Io(«r)  -  2/./  “  IJar)  ]  + 


Cj  [  (2g  ci^  +  — "  )  K^(ar)  +  2p  ^  K^{ar)  ] 
c 

1 

D^2g  i  ^  i  ft  l^(ftt)  -  \  l^iftr)  ]  +  D^2p  i  ^  [  ft  K^{ftr)  + 
j  K^(^r)  ]  ,  (2.9) 


o'^^(r,?,p)  =  2p  [-C^i  ?  «  I^(ar)  +  i  ?  a  K^(ar)  + 

2  2 
,  PP  , 

^  )  l^iftr)  .  D^«  .  -2-  )  K^(/7r)  ]. 

Using  the  notations  introduced  in  (1.36),  the  boundary 
conditions  (2.3)  together  with  (2.9)  the  following  system  of 
linear  algebraic  equations  for  unknowns  C^,  is  found 

C^d^(a)  +  C^d^(a)  -  D^i  d^(a)  -  D^i  d^(a)  =  R^(?,p)  , 

C^d^(b)  +  C^d^(b)  -  D^i  dg(b)  -  D^i  d_^(b)  =  R^(?,P)  , 

-C^i  e^(a)  -  C^i  e^(a)  +  6^(3)  +  e_^(a)  =  0  ,  (2.10) 

-C^i  e^(b)  -  C^i  e^(b)  +  e^Cb)  +  e^(b)  =  0  , 

where 
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dz  dt  = 


00  00  X.V 

R  {?',p)  =  -T  S  H(  c  t  -  |2|  )  e 

-00  0  1 

2  X  V 


-i?z  g-pt 


2  _  2  2 
P  +  ? 


(  2 . 1 : 


The  solution  of  equations  (2.10)  can  be  written  in  the  form 

(2.12) 


T 

=  ¥ 

^^2 

'  ^2  =  ^ 

11 

Q 

T 

II 

where 

( a) 

d,  (a) 

-id^(a) 

-id_^{a) 

djb) 

d2(b) 

-id3(b) 

•id_^(b) 

T  =  det 

-ie^ (a) 

-ie^ (a) 

63  (a) 

e^{a) 

■■ 

-ie^ (b) 

-ie^(b) 

e3(b) 

ejb) 

R^(?,P) 

d2(a) 

-id3(a) 

-  i  d  (  a ) 

R^(?,P) 

d,(b) 

-id3{b) 

-id^(b) 

♦ 

=  det 

0 

-ie^ (a) 

e^(a) 

♦ 

0 

-ie^(b) 

e_^{b) 

d^  (a) 

R^(?,P) 

-id3(a) 

-  i  d_^  (  a ) 

d^  (b) 

R^(? ,P) 

-id3(b) 

~id^(b) 

=  det 

-ie^ (a) 

0 

e  (a) 

e^(a) 

J 

-ie^ (b) 

0 

e  (b) 

4 

' 

d^  (a) 

dj  (a) 

R^(?,P) 

-id_^{a) 

d^(b) 

d^Cb) 

R^(? ,P) 

-id^{b) 

T3  =  det 

-ie^ (a) 

-ie^ (a) 

0 

e^(a) 

. 

.  -ie^(b) 

-ie^ (b) 

0 

e^(b) 

(2.13) 
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d^(a)  -idgla)  R^(?.P) 
d^(b)  -id^(b)  R^(?,p) 

-ie^(a)  e^(a)  0 

-ie^(b)  e^(b)  0 

and  d^  (  .  )  ,  e^(.)  /  '•  =  1  .2  ,3  ,4  ,  are  defined  by  equation 

(  1.36)  . 

Substituting  (2.12)  into  (2.7)  and  using  (1.7)  and  (1.12)  after 
some  manipulations  one  obtains 

— F  ^ 

u;(r,?,p)  =  -  [  T^cx  I^(ar)  -  c  K^(ar)  -  T^i  ?  I^((1r)  - 

i  ?  K  (/?r)  ]  ,  (2.14) 

4  X 

“F  ^ 

u^(r,?,p)  =  ^-  [  i?  I^(c<r)  +  i?  K^(«r)  + 

T3  n  l^inr)  -  /?  K^(^r)  ]  , 

and 

— F 

,P)  =  t  "'^1  i  ?  «  («r)  +  i  ?  a  (or)  + 

2  2 
,  6>p  , 

T3  (?  +  ^  "  2a:  )  ^  ' 
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T  Xp 

o'_^(r,?,p)  =  [  {2p  +  --  )  Ijj(ar)  -  2#^  ?  I^(c<rj 

’  ^  c 

1 

T  Xp^ 

[  (2p  +  ---  )  K^(«r)  +  2p  -  KJar)  ]  - 

c 

1 


2p  i  ?  [  /?  I^(/5r)  -  i  I^(^r)  ]  + 


2^'  i  ?  [  f3  K^(/5r)  +  i  KJ^r)  ]  , 


(2.15) 


0’^^(r,?,p)=  -  [  (pp^  -  2p  a^)  I^(«r)  +  (pp^  - 

2p  o.^)  K^(o(r)  +  Tg  2p  ?  i  r?  I^(/?r)  -  2p  ?  1  K^(r?r)]  . 

Equations  (2.14)  and  (2.15)  represent  the  solution  , called  the 
formal  one  ,of  the  considered  impact  problem  expressed  in  terms 
of  Fourier  and  Laplace  transforms.  The  inversion  of  these 
transforms  would  be  difficult  to  carry  out  because  T,  , 
t.  =1,  2,  3,  4,  are  the  functions  of  transform  parameters  f  and 
p,  equations  (2.13),  (1.36),  (1.28).  Using  the  inversion 

theorems  of  the  Laplace  and  Fourier  exponential  transforms  the 
solution  can  be  written  in  the  form 

1  00  ?-  +  ioo  _  iez  Pt 

u^(r,z,t)  = - -  f  f  u^^(r,?,p)  e  e^  dp  d?"  , 

4  n  i  -oo  ;^-ioD 


00  y+ 100 


P  ^pt 


u  (r,z,t)  = - /  f  u  (r.?,p)  e  e^  dp  d^  , 

Z  .  2  .  .  Z 

4  n  1  -00  ^-100 


00  lOD 


(r,  z,  t)  = 


- ^  /  /  c^;^(r,?,p)  e"''^'  e^"  dp  d? 

4  »T  i  -oo  j^-ioo 


00  y  +  ICO 


o'__  (  r ,  z ,  t )  = 


z  ,,Pt 


/  /  (r,?.p)  e  ^  e^"  dp  d?  , 


4  rr  i  -oo  /'-ioo 
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1 


(  r, z, t ) 


2  . 
n  1 


00  ^  +  ICO 

/  / 

-00  ;^-ioo 


(  r , ,  p )  e 


Uz 


dp  df 
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